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Terahertz (THz) technology has emerged for biomedical applications such as scanning, molecular 
spectroscopy, and medical imaging. Although a thorough assessment to predict potential concerns has to 
precede before practical utilization of THz source, the biological effect of THz radiation is not yet fully 
understood with scant related investigations. Here, we applied a femtosecond-terahertz (fs-THz) pulse to 
mouse skin to evaluate non-thermal effects of THz radiation. Analysis of the genome-wide expression 
profile in fs-THz-irradiated skin indicated that wound responses were predominantly mediated by 
transforming growth factor-beta (TGF-P) signaling pathways. We validated NFkBI- and Smad3/ 
4-mediated transcriptional activation in fs-THz-irradiated skin by chromatin immunoprecipitation assay. 
Repeated fs-THz radiation delayed the closure of mouse skin punch wounds due to up-regulation of TGF-p. 
These findings suggest that fs-THz radiation initiate a wound-like signal in skin with increased expression of 
TGF-P and activation of its downstream target genes, which perturbs the wound healing process in vivo. 

Terahertz (THz) radiation consists of electromagnetic waves that propagate at frequencies from 0.1 to 10 X 
10 12 Hz, the THz range. This radiation is present in our environment as a part of the solar ray spectrum 1,2 . 
THz radiation shows both optical particle-like and electric wave-like characteristics, which has opened new 
fields of research in physics, chemistry, and medicine 2-4 . THz techniques are now used for the development of 
clinical imaging 5 , cancer diagnostics 6,7 , security systems 8 and telecommunication". As THz wave is widely exposed 
to human body with a variety of purposes, we need to characterize the biological response against THz radiation 
in vivo. 

The European project THz- BRIDGE investigated the potential damage of electromagnetic radiation in the THz 
spectral range to biological systems (http://www.frascati.enea.it/THz-BRIDGE/reports/THz-BRIDGE%20Final% 
20Report.pdf). THz radiation in the range of 1 to 3 THz, at 0.45 J/cm 2 , failed to alter the cellular activity or 
differentiation of human keratinocytes, compared to un-exposed cells 10 . Two independent reports have similarly 
failed to identify genotoxic effects of THz radiation in human peripheral blood lymphocytes 11-12 . Recent studies 
using in vitro systems have shown that THz radiation has non-thermally induced impacts on the DNA stability 13-16 , 
which would cause chromosomal aberrations in human lymphocytes 17 or alterations in gene expression with 
accelerated differentiation of mouse stem cells 14-16 . In particular, Titova et al. applied artificial human 3D skin 
tissue model and exposed samples in vitro to broadband THz with pulse energy up to 1 (J.J to detect the signs of 
DNA damage in THz exposed artificial skin tissue 18 . 

In this study, we undertook an integrated bioinformatic and functional analysis to identify genetic alterations 
and following reactions in vivo by THz radiation (Fig. 1 A). Unsupervised approach using mRNA microarray was 
applied to screen THz-responsive genes compared to sham exposed samples. Comparative analysis of the 
expression profiles showed that THz radiation was mostly similar to wound stimulus, not to burning, neutron 
irradiation or UV exposure. This confirmed the in vitro model with artificial skin tissues 18 . Further in silico 
analysis with the differentially expressed genes (DEGs) provided molecular signature responsive to THz irra- 
diation and we found that the wound healing associated signal was predominantly activated via NFkBI- and 
Smad3/4-mediated TGF-P signaling pathway. To verify such a mechanism, we exposed THz repeatedly on 
wounds using an in vivo wound model. Interestingly, we found that over-activated TGF-P signaling with the 
hyper- inflammatory response delayed the healing process of wounds in THz- irradiated mouse skin. 
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Figure 1 | Functional characteristics of the in vivo response to fs-THz radiation. (A) A schematic of the procedures for in vivo exposure and analysis of 
the effects of fs-THz radiation. (B) Differentially expressed genes (DEGs) in THz radiation-exposed skin. Green and red squares with blue lines 
denote acceptable filtering criteria (FCd researcP<0.05). See Table SI for the list of DEGs. (C) Top 10 enriched biological functions associated with 
the 149 DEGs. Each enrichment score is represented by a horizontal orange bar. The vertical blue line indicates acceptable significance (P = 0.05). 
(D) Meta-analysis of the expression of 149 DEGs, compared against gene expression in mouse skin exposed to a variety of stimuli, including UV exposure 
(GSE15618), burn (GSE460), neutron irradiation (GSE25343), and wound (GSE23006). See Fig. S3 and Methods for more detailed information. 



Results 

fs-THz radiation does not affect expression of Hsp70 or histology 
of exposed mouse skin. C57BL/6J mice were exposed to femtosec- 
ond (fs)-THz radiation with a pulse duration of approximately 310 fs 
[full width, at half maximum (FWHM)] and energy of approximately 
0.26 nj/pulse (Fig. S1A and SIB). The frequency spectrum, using 
Fourier transform, ranged up to 2.5 THz (Fig. SIC), at an average 
power density of 0.32 u.W/cm 2 . The accumulated pulse energy for 1 
hour was up to 1.15 mj/cm 2 . Within the measurement error of the 
device (± 0.05°C), there was no change in temperature of the skin of 
C57BL/6J mice that were exposed to fs-pulsed THz radiation 
(Fig. S2A). To evaluate for the presence of THz radiation-induced 



thermal or nonspecific stress, we measured expression of heat shock 
protein 70 (Hsp70) 19 ' 20 at the mRNA (Fig. S2B) and protein levels 19 ' 20 
(Fig. S2C). Thermal stimulus on NIH-3T3 fibroblasts was check as a 
positive control with Hsp70 members including Hspala, Hspa4l and 
Hsphl. We did not observe any notable changes in expression of 
Hsp70 or in the histology of THz-irradiated versus sham-irradiated 
skin of C57BL/6J mice (Fig. S2D). These results indicate that we 
could mine further the non-thermally induced biological consequ- 
ences by THz radiation with the adopted exposure system. 

Characterization of the molecular responses to fs-THz radiation. 

We used microarrays to compare the gene expression profile of 
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mouse skin 24 hoiurs after exposure to sham or fs-THz radiation. 
Through a bioinformatic analysis, we identified 149 differentially 
expressed genes (DEGs) with a mean fold change of signal 
intensity s 1.5 (t-test, P<0.05). Within this gene list were 82 up- 
regulated and 67 down-regulated genes (Fig. IB; Table SI). Gene set 
enrichment analysis with the Ingenuity Pathway Analysis software 
identified the top 10 bio-functions with the significance index of 
Fisher's exact test in THz- irradiated mouse skin (Fig. 1C and Table 
S2). The most significant biological response of 24 hours after THz 
radiation was analyzed as the "healing" process. 

Next, we compared expression of the 149 DEGs in fs-THZ-irra- 
diated mouse skin to gene expression, in datasets from Gene 
Expression Omnibus (GEO), in the skin of mice following exposure 
to various physicochemical stimuli, including after a skin wound 
(Fig. ID and Fig. S3). Normalized datasets were hierarchically clus- 
tered with single-linkage in Euclidean distance, and each of Pearson's 
correlation coefficient (Pec) corresponding to the THz dataset was 
calculated with a two-tailed test of significance (P<0.05). The 
expression of the 149 genes in fs-THz-irradiated skin was mostly 
similar in wound-induced mouse skin, but not in skin exposed to 
other stimuli, such as burning, ultraviolet (UV) exposure, or neutron 
irradiation. 

THz radiation induces expression of genes involved in wound 
response. To validate our results from the analysis of gene 



expression microarrays, we assayed expression of 7 genes from the 
top two functional categories, wound healing (Bmp2, Cd44, Krt6a, 
Lep, Serpinel, Sprrbl, and Thbsl) (indicated in Fig. ID and Fig. S4D) 
and tissue development (Ccll9, Hmgbl, Myf6, Nfatcl, Nr4al, Nr4a3, 
and Vdr) (Fig. 1C and Table S2), by real-time RT-PCR at 1, 12, 24, 
and 36 hours after fs-THz radiation in two mouse strains, C57BL/6J 
(Fig. 2A), BALB/c nude (Fig. 2B and Fig. S4). The expression of 
selected wound response genes were also changed in punch wound 
skin of C57BL/6J mice (Fig. 2C). The expression pattern was largely 
consistent with microarray data (Fig. S4D) and expression of most of 
the wound response genes was increased by 1 hour after exposure to 
THz radiation. Increased expression of Bmp2, Cd44, and Thbsl 
proteins after fs-THz radiation was verified by immunohistoch- 
emical staining (Fig. 2D). Our results establish that the expression 
pattern induced by THz radiation is analogous to that by wound 
stimulus. 

Activation of TGF-P signaling in fs-THz-exposed mouse skin. 

Since the transcription of wound response genes is regulated in 
part by cytokines such as TGF-(3 2123 , we verified activation of 
TGF-P signaling in THz radiation-induced wound responses. First, 
we detected the enrichment of TGF-P signaling pathway in 149 
DEGs by IPA analysis (Fig. 3A). Although TGF-P signaling was 
slightly increased by gene expression microarray, expression of 
Tgfbl mRNA increased at 1 hour after THz radiation, by real-time 
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Figure 2 | Molecular characteristics of responses to fs-THz radiation in mouse skin. (A and B) Relative expression of 7 genes (Bmp2, Cd44, Krt6a, Lep, 
Serpinel, Sprrlb, and Thbsl) by real-time RT-PCR, selected from the top enriched bio-function, "healing" (wound response genes) (See Fig. 1C and Table 
S2), in C57BL/6) (A) and BALB/c nude mice (B). Skin of mice was exposed either to sham or THz radiation, (n = 8 in each group). (C) Expression of 
wound response genes in an in vivo wound model, 24 hours (h) after THz radiation (each, n = 4). (D) Immunohistochemical staining for Bmp2, 
Cd44, and Thbsl. The original magnification used for all images was X 100. A magnified region of staining is shown as an inset in the lower right. 
(A-C, Mean ± standard deviation (SD). *, P<0.05; **, P<0.01; ***, P<0.001). 
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RT-PCR, and decreased thereafter (Fig. 3B). This result was con- 
firmed in BALB/c nude mice and in an in vivo wound model 
(Fig. 3B). As a positive control, we treated NIH-3T3 mouse 
fibroblasts with activators of TGF-P signaling, Activin or TGF-p. 
Similar to our results in THz-irradiated mouse skin, treatment 
with Activin or TGF-P increased expression of Tgfbl and wound 
response genes (Fig. 3C). 

To confirm activation of TGF-P signaling in THz radiation- 
exposed skin, we investigated the activity of TGF-P signaling-related 
transcription factors (TFs). Based on the expression of DEGs in our 
dataset, we identified TFs that could potentially target wound res- 
ponse genes using Ingenuity Knowledge Base software (Fig. 3D). 
Considering most highly activated state (z-score) and its significance 
(P value), we selected Smad3 and NFkBI, that are main downstream 
mediators of TGF-P signaling 24,25 , to examine wound- induced tran- 
scriptional activity in fs-THz-irradiated skin. Next, promoter-wide 
analysis for the transcription factor binding sites (TFBSs) of the 
selected TFs targeting wound response genes (Smad3 for Bmp2, 
Serpinel and Thbsl; NFkBI for Bmp2, Cd44 and Serpinel) was 



carried out using TRANSFAC® professional (Fig. 3E). From in silico 
prediction analyses, we conducted chromatin immunoprecipitation 
(ChIP) assays with primer sets designed to distinguish the promoter 
of each of the genes listed above (Table S4). In mouse skin exposed to 
THz radiation, Smad3 bound to the promoters of Bmp2, Serpinel, 
and Thbsl; NFkBI bound to the promoters of Bmp2, Cd44, and 
Serpinel (Fig. 3E). Smad3 and NFKBlalso bound to the promoter 
of target genes in the skin of THz radiation-exposed BALB/c nude 
mice, wound- induced C57BL/6J mice, and in TGF-P-treated NIH- 
3T3 mouse fibroblasts (Fig. 3E). Although there was not a large 
fraction of promoters bound by the Smad3 and NFkBI TFs, DNA 
binding was quite specific in mouse skin exposed to THz radiation. 

Effect of fs-THz radiation on wound healing. Although we were 
able to identify activation of TGF-P signaling in fs-THz-irradiated 
mouse skin, the irradiated area did not show any histologic evidence 
of wound damage (Fig. S2D). We reasoned that if TGF-P signaling 
was induced in THz radiation-exposed skin without apparent tissue 
damage, it may be affecting wound healing. To this end, we made a 
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Figure 3 | Induction of TGF-P and transcriptional control of wound response. (A) Enrichment analysis of DEGs in signaling pathways. The vertical blue 
line indicates acceptable significance (P = 0.05). "TGF-P Signaling" was analyzed as the most dominant pathway followed by pathways relating to 
translation initiation. (B) Expression of Tgfbl, by real-time RT-PCR, in skin from THz-irradiated C57BL/6J, BALB/c nude mice, and in an in vivo wound 
model. (C) Expression of wound response genes by real-time RT-PCR, in TGF-P-treated NIH 3T3 cells (each, n a 3). (D) in silico mapping of the 
transcriptional regulation for DEGs. Analyzed TFs targeting mRNAs out of database (alphabetized clockwise in outer circle) are placed clockwise from top 
(12:00) in inner circle (See Table S3 for more detailed information). Smad3 and NFkBI are highlighted. (E) Identification for TFBSs of the targeted genes 
regulated by Smad3 and NFkBI . Expression of target genes measured by ChlP-qPCR, in skin from THz-irradiated C57BL/6I and BALB/c nude mice, from 
an in vivo wound model, and in TGF-P-treated NIH3T3 cells. (B, C and E, Mean ± SD. *, P<0.05; **, P<0.01; ***, P<0.001). 
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4-mm punch wound in the skin of C57BL/6J mice and exposed this 
wounded skin to daily fs-THz radiation. After 10 days, punch 
wounds were completely closed in sham-exposed skin. In contrast, 
fs-THz radiation significantly delayed the time for wound closure 
(Fig. 4 A and 4B). The difference in open wound area was also evident 
by analysis of tissue histology (Fig. 4A). We observed slower re- 
epithelialization with bigger scab formation in mouse skin that was 
continuously exposed to THz radiation. The results of a multiplex 
enzyme-linked immunosorbent assay (ELISA) for expression of 
TGF-P family members showed significant elevation of Tgfbl 
protein both at 5 and 9 days following initiation of THz radiation 
to wound-exposed mouse skin (Fig. 4C). Thus, elevation of Tgfbl 
expression may be one explanation for the inhibited wound healing 
process in THz- irradiated mouse skin. These results are consistent 
with previous findings that elevated plasma levels of TGF-P inhibit 
re-epithelialization with scarring 26,27 (Fig. 4D). 

Discussion 

The wide variety of applications for THz technology has raised issues 
of biological safety due to exposure to THz-associated radiation 1 . In 
this study, mRNA expression profiling detected wound responses 
induced by fs-THz radiation in mouse skin. Although we could 
not find either macroscopic or microscopic evidence of tissue 
damage in irradiated skin, our results suggest that THz radiation 
induces wound-like stimuli in keratinocytes through activation of 
TGF-P signaling. 

Wound response is a dynamic process that includes inflamma- 
tion, granulation, re-epithelialization, matrix formation, and tissue 
remodeling 28 . Keratinocytes have a role in re-epithelialization as a 
consequence of proliferation and differentiation in basal layers of the 



epidermis. Although the process of early inflammation is mediated 
by fibroblasts and macrophages, keratinocytes up-regulate express- 
ion of cytokines to recruit neutrophils to sites of injury 22,23 . Among 
the up-regulated cytokines, we discovered significant up-regulation 
of TGF-P signaling 24 hours after mouse skin was exposed to THz 
radiation. TGF-P signaling has pleiotropic effects on inflammation 2 *. 
Others have shown that elevated levels of TGF-P in keratinocytes 
significantly inhibit wound re-epithelialization in transgenic mice 
overexpressing TGF-p 27 . We have similarly shown that repetitive 
exposure to THz radiation results in elevated levels of TGF-P and 
delayed closure of punch wounds in mouse skin. 

TGF-P signaling is primarily propagated by Smads and a complex 
of NFkB 21,24,25,30 . We have shown that THz-sensitive genes are mostly 
associated with wound response and tissue development. We have 
further shown that transcriptional activation of these genes is 
mediated by the transcription factors Smad3/4 and NFkBI. In 
Fig. 4D we propose a mechanism by which THz radiation promotes 
a wound response in skin. However, the initial physical interactions 
between THz radiation and skin that trigger this signaling cascade 
remain unclear. By abrogating TGF-P-induced stimulation of col- 
lagen gene expression, myofibroblast transdifferentiation, and 
Smad-dependent activity in normal fibroblasts, PPAR-y may play 
a physiologic role in regulating the pro-fibrotic response 31 . We have 
shown increased levels of cytokines, such as TGF-P, following expo- 
sure of skin to THz radiation. 

We found that a number of genes related to would response, such 
as Bmp2, Cd44, Thbsl, and Serpinel, were specifically altered by THz 
radiation. UV radiation has also been reported to increase expression 
of TSP-1 at the mRNA and protein levels, in human skin dermis 
and dermal fibroblasts 32 . The transmembrane glycoprotein CD44 
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Figure 4 | Delayed wound healing associated with fs-THz radiation. (A) Gross and microscopic time course photographs of skin wounds following 
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model of wound response induced by fs-THz radiation. TGF-P signaling is activated in THz-exposed skin, which triggers transcription of wound healing- 
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is currently thought to be the main cell surface receptor for the 
glycosaminoglycan hyaluronate. Although UV radiation reduced 
expression of CD44 and hyaluronate in the epidermis of hairless 
mice, topical retinol or retinoic acid increased the basal level of 
epidermal CD44 33 . The responses to a pulse of THz radiation may 
overlap with responses to a series of environmental stresses such as 
UV, LED, and retinol. 

There have been a few reports suggesting direct interactions 
between THz radiation and biomolecules, which regulate the effects 
of THz radiation. The structure and function of biomolecules are 
dynamically controlled via their interactions with water 34 ' 35 . Based on 
the fundamental principle that the stretching and bending of hydro- 
gen bonds is affected in the THz domain 36 , it is conceivable that the 
vibrational, torsional, and librational dynamics of molecules near 
water are affected by exposure to THz radiation. Alexandrov et al. 
also proposed that THz radiation may significantly perturb the nat- 
ural state of dsDNA 13 . Thus, THz radiation may alter important 
molecular processes involved in gene expression and DNA replica- 
tion. Others have demonstrated augmented expression of adipogenic 
differentiation genes (Adiponectin, Glut4, Fabp4, and Pparg) in 
mouse stem cells after 9 hours of THz exposure 14,15 . In our experi- 
mental conditions, within 1 hour of exposure of mouse skin in vivo to 
THZ radiation, prediction analysis of putative transcription factors 
governing DEGs showed slight activation of PPAR-y (Table S3). 
Unfortunately, although we noted a transient change in the express- 
ion of Adiponectin, Glut4 y Fabp4, and Pparg 12 hours after THz 
radiation (data not shown), these four genes were not detectable after 
filtering for DEGs at 24 hours after THz radiation. 

In summary, we demonstrated for the first time that THz radiation 
produces cellular responses against wound-like stimulation. Although 
THz radiation barely penetrated the outermost region of the epidermis, 
absorption of THz radiation increases with frequencies in the THz 
domain 37 . Overall, our results suggest that accumulated THz exposure 
may alter hydration in skin sufficiently to alter molecular signaling 
cascades. We have observed non-thermally induced effects of fs-THz 
radiation in mouse skin. Considering the rapidly increasing application 
of THz in biology and biomedicine, it will be increasingly more import- 
ant to understand the potential health hazards of THz radiation. To 
this end, our results are particularly important as they may be used as 
biomarkers for the detection of THz radiation exposure. 

Methods 

In vivo fs-THz exposure. We used a fs-THz exposure system which delivers optical 
pulses produced in a 3.5 W Ti:Sapphire regenerative laser amplifier with a center 
wavelength of 800 nm and a repetition rate of 1 kHz 3839 (Fig. SI A). Briefly, the pulse 
width was approximately 180 fs (FWHM). This was delivered as approximately 1 mj 
of energy out of the regenerative amplifier. The 800 nm laser beam was collimated 
and split into two parts by a beam splitter. One part {up to 750 uj) was delivered to a 
10 X 10 mm ZnTe crystal with a thickness of 1 mm to generate fs-THz pulses via 
optical rectification. This THz electric filed was mapped out by electro-optic sampling 
method in another 1 m thick ZnTe nonlinear crystal. The other part of the 800 nm 
laser beam with a pulse energy (up to 30 uj) was used as a gate pulse. This gate pulsed 
was delayed using an optical delay line and collinearly focused with the THz beam 
onto the ZnTe crystal. In electro-optic sampling, the THz beam induces birefringence 
in an electro-optic medium. If the polarized gate pulse traveled through the ZnTe 
crystal at the same time with THz pulse, the induced birefringence rotated the 
polarization of the gate pulse. The magnitude of rotation was proportional to the THz 
electric field strength. The signal was collected by a combination of a quarter-wave 
plate, a Wollaston polarizer, and a pair of balanced photodiodes. The signal current is 
then sent to a lock-in amplifier {SR830; Stanford Research System, Sunnyvale, CA). 

The entire THz waveform was measured by gate pulses as a function of optical 
delay time {Fig. SIB). The amplitude spectrum was obtained after Fourier trans- 
formation of the THz field, which covered frequencies up to 2.5 THz (Fig. SIC). 
Laser-based, THz pulse energy (approximately 0.26 nj/pulse) was measured with a 
Golay cell (Microtech Instrument, Inc., Eugene, OR). In order to minimize the 
absorption by water vapor, a region of THz beam path was enclosed and purged with 
dry air. A THz low pass filter was used to remove the extra visible source while this 
experiment was in progress. In the THz-targeted area, the accumulated pulse energy 
over 1 hour was up to a maximum of approximately 1.15 mW/cm 2 . A collimated THz 
beam of energy was delivered with a 800-nm laser, targeted to a 1 0-mm-diameter area 
of skin on the dorsum of a mouse that was secured within an acrylic mouse restrainer. 
Temperature in the target area of sham- and THz-exposed skin was monitored with a 



high -resolution infrared camera (T335; FLIR Systems, Boston, MA). Real-time 
temperature data was analyzed using FLIR QuickPlot software (FLIR Systems, Fig. 
S2A). 

Animals and wound model. Six- to eight- week- old male C57BL/6J and BALB/c nude 
mice were used in this study. All in vivo experimental procedures, as well as use and 
care of animals was in compliance with requirements of the Seoul National University 
IACUC (approval number: SNU- 101108-1, SNU- 120112-3). Temperature (22 ± 
2 LI C) and humidity {50 ± 10%) were automatically maintained with an air controlling 
system in mouse housing. Mice were anesthetized by intraperitoneal injection of a 
cocktail of Zoletil (30 (ig/g) and Rompun (10 ug/g) just before exposure to THz 
radiation. A marked area of the shaved dorsal skin of mice was used to administer 
THz radiation. After a 1 hour exposure to THz radiation, a full thickness of exposed 
skin was removed, at indicated times, for analysis. 

Four pairs of wounds were made on the shaved dorsal skin with a 4-mm sterile 
biopsy punch (Kai Medical, Honolulu, HI). For the experiments in which THz 
radiation was delivered to wounded mice, we spread depilatory cream on the shaved 
back of mice. Images of each wound were acquired by a digital camera before radi- 
ation with sham- or THz radiation for 1 hour daily. Two mice in each of two groups 
were sacrificed after obtaining an 8 mm punch biopsy at 1, 3, 5, 9, and 12 days 
following radiation. To prevent experimental error from mice rubbing their wounded 
skin against each other, each test and control-treated mouse was housed in its own 
cage. In addition, bedding, food, and water for mice were changed every other day. 

Size-compensated photographs of wounds at indicated time points were examined 
manually to determine wound size and converted to binary images with ImageJ 
software (http://rsb.info.nih.gov/ij). The percentage of closed wound area for each 
mouse was averaged, and standard deviation was calculated and displayed in figures 
as error bars. The difference in the size of the wounds between sham- and THz 
radiation -exposed groups was determined by a two-tailed Student's l-test (P<0.05). 

mRNA microarray and bioinformatic analysis. Total RNA extracted from mouse 
skin was labeled and hybridized on the Mouse Gene 1.0 ST Array (Affymetrix, Santa 
Clara, CA), according to the manufacturer's protocol. The expression profiling have 
been deposited at the Gene Expression Omnibus (GEO) with the accession number 
GSE 44671. THz-responsive genes were filtered out of global gene expression using 
moderated i-statistics, based on an empirical Bayesian approach and a Benjamini- 
Hochberg false discovery rate (BH FDR)-adjusted P value. Statistical criteria to 
identify the differentially expressed genes (DEGs) were defined by an absolute fold- 
change (^1.5) and P value (<0.05). Statistically significant DEGs were examined for 
possible interactions with other molecules, canonical signaling pathways, and bio- 
functions with the Ingenuity Pathway Analysis software package {IPA; Ingenuity® 
Systems), using standard protocols 40 ' 41 . 

To compare DEGs from THz-irradiated skin with DEGs identified when mouse 
skin was subjected to a variety of exposures, we obtained the following gene 
expression datasets from NCBI GEO: GSE15618 {UV exposure), GSE460 (burning), 
GSE25343 (neutron irradiation), and GSE23006 (wounding). Data was normalized 
according to each of authors' methods. The fold-change in mRNA expression was 
calculated using the difference of the average values to each control. Hierarchical 
clustering trees were generated with single -linkage in Euclidean distance, using the 
MeV (MultiExperiment Viewer, http://www.tm4.org). A two-tailed Pearson's cor- 
relation was applied, with the threshold for significance set at P<0.05. 

Real-time reverse transcriptase-polymerase chain reaction (RT-PCR). To 

quantitatively estimate the mRNA expression of target genes, cDNA was synthesized 
using a ReverTra Ace® qPCR RT kit (Toyobo, Osaka, Japan). PCR was performed in a 
Mx3005P QPCR System (Agilent Technologies, Santa Clara, CA) using SYBR® 
Premix Ex Taq 1M (Takara Bio Inc., Shiga, Japan), according to the manufacturer's 
instructions. The mean SCt (change in cycle threshold, Ct) was normalized to the 
mean Ct for the housekeeping gene, Gapdh. PCR reactions were performed in 
triplicate for each sample; error bars show standard deviation. Statistical significance 
was determined using a Student's f-test. P values are two-tailed and significance was 
accepted at P<0.05. The sequence information for primer sets to validate mRNA 
expression is available in Table S4. 

Immunohistochemistry. Biopsied skin tissues were deparaffmized and stained with 
hematoxylin and eosin, or hybridized with primary antibodies corresponding to each 
secondary antibody. Antibodies for Bmp2, Cd44, Thbsl, and Hspala were purchased 
from Santa Cruz Biotechnology Inc. {SANTA CRUZ, CA). Stained images were 
acquired using a microscope (DM2500; Leica Microsystems, Wetzlar, Germany) 
mounted with a CCD camera {DFC290; Leica Microsystems). Image analysis was 
performed with Leica Application Suite (Leica Microsystems) software. 

Cell culture and TGF-p treatment. NIH 3T3 mouse fibroblasts were obtained from 
AATC (Manassas, VA), and maintained as the protocol provided by AATC. Heat 
shock at 38.5, 40 and 43 °C were given for 1 hour and placed to 37 IJ C incubator. For 
TGF-p treatment, we first treated cells with 10 uM SB431542 (Tocris, Bristol, UK) for 
24 hours, and added 10 ng/ml Activin (R&D Systems, Minneapolis, MN) for 2 hours, 
or 2 ng/ml TGF-J5 (R&D Systems) for 12 hours. 

In silico analysis of transcription factors. Upstream TFs regulating DEGs identified 
from our analysis of microarray data were predicted by core- analysis using the IPA 
software package. Among TFs identified by the IPA software were NFkBI and 
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Smad3/4. These TFs were studied further because of their high z-scores and 
association with TGF-|3 signaling in wound -inducible responses 21,24 ' 25 . Interactions of 
these two TFs with their target mRNAs were highlighted using Path Designer, in the 
IPA software package. 

Chromatin-immunoprecipitation (ChIP). Transcriptionally-regulated DNA sites 
for target genes of NFkBI and Smad3/4 were analyzed with TRANSFAC® 
professional (BIOBASE, Wolfenbuttel, Germany), following standard protocols 42 . 
Promoter sequences for each target gene were identified with Regulatory Sequence 
Analysis Tools (RSAT; http://rsat.ulb.ac.be), with an Ensembl ID. Thousands of 
promoter sequences were retrieved with the TRANSFAC® database. Potential TF 
binding sites in the promoter of target genes were scanned. Two candidate sites were 
selected using the highest matrix similarity and core similarity scores. 

Skin tissue sections and cells were fixed in 1% formaldehyde at room temperature. 
Genomic DNA was extracted from tissue or cells using a Pierce® Agarose ChIP kit 
(Thermo Scientific, Waltham, MA). Instead of sonication, we applied micrococcal 
nuclease digestion to fragment DNA into a uniform size. ChlP-grade anti-NFKBl 
(Abeam, Cambridge, UK) and anti-Smad3 (Abeam) antibodies were used for immu- 
noprecipitation (IP) with the digested chromatin. Finally, purified DNA samples, 
including input and IP, were prepared for qPCR detection with the specific primers 
listed in Table S4. IP with negative control IgG or positive control RNA polymerase II 
antibody and all the designed promoter primers, as well as primers for Gapdh, were 
validated by PCR amplification and agarose gel electrophoresis. Real-time PCR was 
carried out with a Mx3005P QPCR System using SYBR® Premix Ex Taq™, according 
to the manufacturer's instructions. PCR reactions were performed in triplicate, nor- 
malized to the ACt value of IP to that of input, and expressed as the percentage of input, 
with error bars representing standard deviation. A two-tailed Student's Mest was used 
to determine statistical significance; significance was accepted at P<0.05. 

Measurement of TGF-p. The protein level of TGF-|3 in mouse skin was measured 
using a Bio-Plex Pro IM Assays kit, according to the manufacturer's protocols (Bio- 
Rad, Hercules, CA). Lysates from mouse skin were adjusted to a uniform 
concentration before measurement in a Bio-Plex® 200 System (Bio-Rad). Beads for 
anti-TGF- P were mixed with lysates, and the fluorescent intensity was converted into 
a concentration of TGF-p. Each value was an average of triplicate experiments with 
standard error. The difference of TGF-P concentration between sham- and THz 
radiation- exposed skin was considered to be significant for P<0.05. 
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